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Abstract 



Algoritlims are described and Maple implementations are provided for finding all quandles of 
order n, as well as computing all homomorphisms between two finite quandles or from a finitely 
presented quandle (e.g., a knot quandle) to a finite quandle, computing the automorphism group 
of a finite quandle, etc. Several of these programs work for arbitrary binary operation tables 
and hence algebraic structures other than quandles. We also include a stand-alone C program 
which finds quandles of order n and provide URLs for files containing the results for n = 6, 7 
and 8. 
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1 Introduction 

In 1980, David Joyce introduced a new algebraic structure dubbed the quandle. Quandles are tailor- 
made for defining invariants of knots since the quandle axioms are essentially the Reidemeister moves 
written in algebra. Associated to any knot diagram, there is a quandle called the knot quandle which 
is a complete invariant of knot type up to homeomorphism of topological pairs. 

The history of quandle theory is a story of rediscovery and reinvention. Quandles and their 
generalization, racks, have been independently invented and studied by numerous authors (H], H, 
[S], ^U], etc.) and classification results for various subcategories of quandles have been obtained 
by various authors (0, In [Z|> the third listed author and a coauthor described a way of 

representing finite quandles as matrices and implemented algorithms for finding all finite quandles, 
removing isomorphic quandles from the list, and computing the automorphism group of each quandle. 
As we later learned, some of our work has duplicated the efforts of others 

This paper is intended to reduce future duplication of effort by describing the algorithms for 
computation with finite quandles implemented in "T* and other recent projects, as well as an improved 
algorithm for finding quandle matrices. The C source for our implementation of this algorithm as well 
as Maple implementations of algorithms for computing with finite quandles and the lists of quandle 
matrices of order 6, 7 and 8 are available for download at http://www.esotericka.org/quandles. 
Additional maple code corresponding to current and future projects will be made available at the 
same site, such as an algorithm for finding all Alexander presentations of a quandle when such exist 



nn. 



2 Quandles, quandle matrices, and homomorphisms 



Definition 1 A quandle is a set Q with a binary operation > 



Q X Q ^ Q satisfying 



(i) for every a; e Q we have x > x = x, 



(ii) for every x,y G Q there is a unique z ^ Q such that x 



z> y, and 



(iii) for every x,y, z £ Q we have (x> y) t> z = (x > z) > (y t> z). 
If (Q,t>) satisfies (ii) and (iii), Q is a rack. 



Axiom (ii) says that t> is right-invertible; for every y £ Q, the map fy ■ Q ^ Q defined by 
fy{x) = xt>y is a bijection (indeed, a quandle automorphism). Denote the inverse map fy^{x) = x<iy. 
Then {Q,<i) is also a quandle, called the dual of {Q,>); not only is < self-distributive, but it is an 
easy exercise to check that t> and <i distribute over each other. 

Standard examples of quandles include groups, which are quandles under conjugation g > h = 
h^^gh as well as n-fold conjugation g>h = h~"'gh'^, denoted Conj(G') and Conj„(G) respectively, 
and Alexander quandles, which are modules over the ring A = Z[t^^] of Laurent polynomials in one 
variable with integer coefficients, with quandle operation given by 



A finite quandle Q may be specified by giving its quandle matrix Mq, which is the matrix 
obtained from the operation table of Q — {xi,X2, . . . ,a;„} (where the entry in row i column j is 
Xi >Xj) by dropping the xs and keeping only the subscripts. In 7. it is noted that, unlike arbitrary 
binary operation tables or indeed even rack tables, quandle axiom (i) permits us to deduce the 
column and row labels from the elements along the diagonal of a quandle matrix. 

Example 1 Let Q = R4, the dihedral quandle of order 4, which has underlying set Q ~ {xi — 
0, X2 — l,X3 — 2, X4 — 3} with quandle operation Xi>Xj — x^j-i (mod4)- Then Q has operation table 



Xi X2 X^ X4 

Xl X\ Xj, X\ X-i 

X2 X4 X2 Xi X2 and hence matrix Mji^ = 

X3 X3 Xl X3 Xl 

X4 X2 X4 X2 X4 



A map (j) : Q Q' from a quandle Q = {xi, . . . , Xn} to a quandle Q' — {yi, . . . , ym} rnay be 
represented by a vector v = {(p^xi), 4>{x2), ■ ■ ■ , (t){x„)) G Q'". Such a vector v then represents a 
homomorphism iff (f>{xi > Xj) = 4>{xi) > 4>{xj), that is, iff we have 



for all Xi,Xj e Q where A — Mq, B — Mqi, and the notation M[i,j] indicates the entry of M in 
row i column j . 

In presentations of quandles by generators and relations are defined. In [Tl', it is observed 
that all finitely presented quandles may be written with a short form presentation in which every 
relation is of the form a = bo c where o G {>, <}. In particular, a knot quandle has a presentation 
with n such short relations where n is the number of crossings in the diagram. Moreover, we may 
assume (rewriting if necessary) that every relation is written in the form a = b> c and that no two 
relations of the form a = bc> c and a' = b> c are present, since if a — b> c and a' = bt> c are both 
present we can replace every instance of a' with a and remove the generator a' without changing 
the presented quandle; in particular, if our quandle is a knot quandle, Reidemeister type I moves^ 
induce such a replacement. 

Definition 2 Let Q = (I, 2, . . . , 71 | ai = &i >ci, . . . , a„i — bm^Cm, m < n^) he a. short form quandle 
presentation such that no two relations of the form a.; ~ bi t> Ci and aj = bi t> Ci with a.; ^ aj are 
present. The matrix MP e A/„(Z) with 



is the matrix of the presentation Q. Note that a quandle matrix for a finite quandle is the matrix 
of a presentation of a finite quandle, so this definition generalizes the notion of quandle matrices to 
finitely presentable quandles. 

^Reidemeister moves are described in 1141 and many other works. 



X > y — tx + {1 — t)y. 



v[A[i,j]]^B[v[i],v[j]] 




k = i> j a listed relation 
otherwise 



Example 2 

3 

2 

1 



The pictured trefoil knot diagram has quandle presentation (1, 2, 3 |1 = 2 > 3, 2 = 3 > 1, 3 = 1 > 2). 
The relations are determined at a crossing by looking in the positive direction of the overcrossing 
strand indicated by the given orientation; the relation is 

(left-hand undercrossing) = (right-hand undercrossing) i> (overcrossing). 

See 0] or for more. 

This matrix representation gives us a convenient way to do computations involving quandles, 
including the quandle counting invariant for knot quandles or other short form quandles with respect 
to a finite target quandle. The next section describes algorithms for doing computations with 
quandles and refers to implementations in Maple T7 and C [S]. 

3 Algorithms 

The goal of the computations in 171 was to find all quandles of a given finite order. Originally, we 
wrote separate programs for each value of n; (T^ includes one example of such an implementation, 
quandleslistS. We later wrote a more general program which works for arbitrary n, though due 
to the large number of columns to be checked, for values of n > 6 we decided to implement a 
stand-alone version suitable for distributed computing. 

The algorithm implemented in qucindleslist takes a number n and generates a list of all n x n 
standard form quandle matrices. A matrix M 6 M„(Z) is a quandle matrix in standard form iff it 
satisfies the following three conditions: 

(i) for i G {1, . . . , n}, M[i, i] = i, 

(ii) every column in M is a permutation of {1, . . . , ti}, and 

(iii) for every triple I < i,j,k < n we have M[M[i,j],k] = M[M[i, k],M[j, k]]. 

To guarantee that conditions (i) and (ii) are satisfied, we start by getting a list of all permutations 
of {1, . . . ,n}. The program listperms takes a number n and produces a list of all permutations 
p £ Yin, represented as vectors [p(l), p(2), . . . , p(n)], in the dictionary order. 

The ith column in a standard form quandle matrix has entry i in the zth position. The pro- 
gram listpermsi takes a pair of positive integers {n, i) and outputs a list of all permutations of 
{1, 2, . . . , n} p E Sn which fix the element i. 

To test quandle axiom (iii), we note that the first time any triple (i,j,k) fails to satisfy the 
axiom, we can exit the program and report that the matrix is not a quandle. This is implemented 
in qStest. 

For a fixed value of n, we can then simply run over a series of nested loops, testing each resulting 
matrix for quandle axiom (iii), since by construction axioms (i) and (ii) are already satisfied. The 
program quandleslistB is an example of this. 

The program quandleslist finds a list of all quandle matrices of a given size n. To find all 
n X n quandle matrices for arbitrary n, qusindleslist finds all control vectors v[i\ with n entries 
using listmaps, a program which takes two inputs a and b and outputs a list of all a-component 
vectors with entries in {1, . . . , 6}. Each entry in the control vector corresponds to a column in the 
output matrix; for each such control vector, an n x n matrix Af [i, j] is produced whose ith column 
is L[n, i][ti[i]], where L[n, z] is the output of listpermsi (n, i) . These matrices are then tested for 




quandle axiom (iii) using qStest. For completeness, we include a program which tests a matrix for 
all three quandle axioms, qtest. 

Since every n-component vector with entries in {!,..., m} can be interpreted as a map from 
{1, . . . , n} to {1, . . . , m}, we can use listmaps to compute the set of all homomorphisms from one 
finite quandle to another. Let A g M„(Z) be an n x n quandle matrix and B € Mm(Z) an m x m 
quandle matrix. Then the vector v £ 17^ , 1 < v[i] < m represents a quandle homomorphism 
v: Bi^ 

v[A[i,i]] = v{i>j) = v{i)i>v{j) = ^[j]], 

as noted in section |21 The program homtest takes two quandle matrices and a vector and reports 
whether the vector represents a quandle homomorphism or not. 

The program homtest handles the case where A is either a finite quandle matrix or a presentation 
matrix for a finitely presented quandle; in the former case, the program simply tests whether the 
assignment of generators {1, . . . , n} in the quandle with presentation matrix A to elements {1, . . . , m} 
in the finite quandle B satisfy the relations defining A by ignoring any zero entries in A. 

We make use of nextmap, a procedure which takes as input a vector v and number n and returns 
the next m-component vector with entries in {1, 2, . . . , m} in the dictionary order, to get a list of 
all homomorphisms from the quandle with matrix A to the quandle with matrix B in the program 
homlist. The program homcount counts the number of homomorphisms from one finite quandle to 
another. If ^ is a knot quandle presentation matrix, then homcount computes the quandle counting 
invariant, i.e., the number of quandle colorings of the knot diagram defining A by the finite quandle 
B. Alternate methods of computing the quandle counting invariant for finite Alexander quandles 
are described in 

After the first version of this paper was completed, we implemented a much faster algorithm 
for finding quandle homomorphisms, homlist2. This program uses a |i3|-component vector with 
entries in {0, 1, . . . , \A\} as a template for a homomorphism, with entries acting as blanks to be 
filled in. The program keeps a working list of such templates, systematically filling in zero values and 
propagating the value through the template using homf ill. The procedure homf ill takes as input 
a quandle matrix B, a quandle presentation matrix A and a template vector v and systematically 
checks every pair of entries for the quandle homomorphism condition j]] = _B[i;[i], v[j]\, filling 

in zeroes where possible and reporting "false" if a contradiction is found. 

Since an isomorphism is a bijective homomorphism, and a bijective map is represented by a 
permutation v : {1, . . . , n} — > {1, . . . , n}, we can test whether two quandles given by matrices are 
isomorphic by running through the list of permutations of order n and testing to see whether any 
are homomorphisms. The program isotest returns "true" the first time it finds an isomorphism 
and "false" if it gets through all nl permutations without finding an isomorphism.^ 

Replacing listmaps in homlist2 with permute (n) and setting B = A gives us the automorphism 
group of the quandle with matrix A, autlist, represented as a list of permutation vectors. 

In [7], a slightly different method of determining the automorphism group of a quandle was used. 
Specifically, permuting the entries of a quandle matrix A by a permutation p applies an isomorphism 
to the defined quandle, but the new matrix now has its rows and columns out of order. To restore 
the order, we conjugate by the matrix of the permutation; the resulting matrix was called p{A) in 
[7]. In particular, a permutation p is an automorphism of A iff p{A) = A. To compute Aut{A) in j7j, 
we ran a loop over the list of permutations given by listperms and noted which ones preserved the 
original matrix A. Here, we include a program stndiso which computes the standard form matrix 
p{A) given a quandle matrix A and a vector v representing the permutation p. 

Finally, once we have a list of quandle matrices of order n, we want to choose a single rep- 
resentative from each isomorphism class. The program reducelist takes a list of quandle ma- 
trices and compares them pairwise with isotest, removing isomorphic copies and outputting a 
trimmed list. The program reducelist works for short lists; an improved algorithm, implemented 
as reducelist2, is better for longer lists, but neither is sufhcient to reduce the rather lengthy lists 
of quandles of order 7 and 8 in a reasonable amount of time. 

faster version of this program using orbit decompositions of finite quandles is described in 



We note that several of these programs, notably homtest, homlist, homlist2, homcount, isotest, 
autlist, and reducelist are not quandle-specific but apply as written to any bmary operation de- 
fined using a matrix as operation table. These facts are exploited in jllj . in which the authors give 
a program which determines all Alexander structures on a quandle, if there are any, using matrices 
to represent the Cayley table of an abelian group. 

We have also implemented a stand-alone version of quandleslist, written in C (see jHI); it 
writes a list of quandle matrices in Maple format to an output file. 

In our initial version of the stand-alone program, several instances of the program could be run 
in parallel on networked machines using a control file to ensure that separate instances do not repeat 
the same computations. However, sufficient improvements were made to the algorithm by pruning 
the search space that the current version can handle the n = 8 case on a single processor, though 
the n = 9 case is still out of reach even with a large network. 

The first improvement was to introduce a partial test versus axiom (iii) after generation of each 
column. In many cases we can find entries that violate the axiom well before the entire matrix is 
generated, which allows vast portions of the search space to be pruned. 

The second improvement was to notice when all of the interior coordinate values as well as the 
left-hand side value of the axiom (iii) equality have been computed, but the right-hand side value 
has not. In this case we can constrain a row of a future column to be equal to the left-hand side 
value. This reduces the number of rows that must be permuted when searching that column, which 
further prunes the search space. The pruning effect is magnified the earlier these constraints are 
added. For example, with n—7 and all else held equal, adding a single constraint to column 3 saves 
(6! — 5!) * ((6!)^) or 2.3 x 10^^ tests, whereas adding a constraint to column 7 saves only 6! — 5! or 
600 tests. 

The effect of the two improvements can be seen in the following tables. It is interesting to note 
that although there is nothing in the program to prevent it (and reasonable amount of code to 
encourage it), we never add constraints to column 2, nor do we ever add more than one constraint 
per column, nor do we detect addition of conflicting constraints. 



n 


5 


6 


7 


8 


total search space 


8.0 X 10« 


3.0 X 10" 


1.0 X 10^" 


4.1 X 10^*^ 


with early testing 


8400 


715680 


1.0 X 10*^ 


n/a 


with forward propagation 


1154 


53500 


5.0 X 10^ 


7.7 X IQS 


total quandles 


404 


6658 


152900 


5225916 



Table 1: Number of complete matrices tested. 



column 


2 


3 


4 


5 


6 


n = 5 


186 


8736 


14626 






n — 6 


4728 


1090404 


8418374 


1187556 




n = 7 


154680 


2.3 X 10^ 


1.8 X IQio 


3.8 X 10^ 


1.8 X 10^ 



Table 2: Number of tests pruned with early testing. 



column 


3 


4 


5 


6 


7 


n — 5 


164 


179 


290 






n ~ 6 


3558 


4396 


3348 


5020 




n = 7 


115872 


228384 


91452 


82910 


117430 



Table 3: Number of times columns constrained beyond axiom (i). 
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